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Higher order and multicomponent generalizations of the nonlinear Schrodinger equation are im- 
portant in various applications, e.g., in optics. One of these equations, the integrable Sasa-Satsuma 
equation, has particularly interesting soliton solutions. Unfortunately the construction of multi- 
soliton solutions to this equation presents difficulties due to its complicated bilinearization. We 
{T") ■ discuss briefly some previous attempts and then give the correct bilinearization based on the interpre- 

tation of the Sasa-Satsuma equation as a reduction of the three-component Kadomtsev-Petviashvili 
hierarchy. In the process we also get bilinearizations and multi-soliton formulae for a two component 
generalization of the Sasa-Satsuma equation (the Yajima-Oikawa-Tasgal-Potasek model), and for a 
(2 + l)-dimensional generalization. 



> 



X 



PACS numbers: 05.45.Yv, 42.65.Tg, 42.81. Dp, 02.30.Jr 



I. INTRODUCTION 



s 

One of the most interesting applications of solitons is in the propagation of short pulses in optical fibers (for an 
overview, see p|). The basic phenomena are described by the nonlinear Schrodinger equation, but as the pulses get 
shorter various additional effects become important. In 2] Kodama and Hasegawa derived the relevant equation with 
higher order correction terms, the generic form of such an equation is (in the optical fiber setting the roles of time 
and space are usually reversed) 

(N , 

>'■ iq i + a 1 q TT + a 2 \q\ 2 q + t[l3 1 q TTT + (3 2 \q\ 2 q T + /3 3 q(\q\ 2 ) T ] = 0, (1) 

in ■ 

where the ai,fti are real constants and q a complex function. The first three terms form the standard nonlinear 
Schrodinger equation (nls) and the /3j terms are the perturbative corrections. Usually one chooses the scaling so that 
| o;2 = 2ol\. In this paper we assume /3i ^ 0. 
(*C) • Our main concern here is the bilinearization and multi-soliton solutions of the Sasa-Satsuma equation (SSnls) 3], 
which is a particularly interesting integrable example in the above class. In this introductory section we discuss some 
basic properties of Q, its integrable special cases and their multicomponent generalizations. In particular we show 
• t-h ' that many previous attempts to solve these equations have produced only rather trivial solutions, in which the complex 
and multicomponent freedom has been "frozen" . The reason for this turns out to be in the incorrect bilinearization 
that was used in those papers. The correct bilinearization (presented in Sec. [n] with detailed derivation in Sec. IIII|) 
follows once we identify SSnls as a reduction of the three-component Kadomtsev-Petviashvili (KP) hierarchy, and 
then we also obtain general multi-soliton solutions. 



A. Gauge transformation 

In order to understand the complex structure of Q it is important to isolate the gauge (phase) invariance and fix 
the gauge. First let us recall that the nls part of Q (i.e., if (3% — 0) is invariant under the combined gauge/Galilei 
transformation 

q(^r) = e i ^- v ^y(x,t), x = r - 2«£, t = £. (2) 
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The full equation JQ) is not invariant under , but if we try the transformation 

g(£, r) = e l( - ClT+C2 ^y(x, t), x = t + c 3 £, t = f , c % real constants, (3) 

we find that if 

C3 = cx{-2a x + 3ftci), c 2 = c?(-ai + cift), (4) 
then Q is /orm invariant: the equation for y(x,t) is as in with ft unchanged, but the oti change according to 

ai — > Si = ai — 3ftci, a2 — » S2 = a2 — ftci. (5) 

We can therefore use this transformation to put a\ = a 2 = 0, provided that 

3ft a 2 = ft>a<i. (6) 

[In the usual normalization ol\ = h, a 2 = 1 © means ft = 6ft.] In all integrable cases (along with some non-integrable 
cases appearing in the literature) JHJ is satisfied, and we assume it from now on. 

On the basis of the above we fix the gauge (J2J by requiring that ai — in the equation, and compare results only 
in that uniquely defined gauge. 



B. Integrable cases 

The integrability of the class of equations has been studied by a number of authors using Painleve analysis 
0- IE 0] and other methods Q , with the consistent result that if ft , ft 7^ there are precisely two integrable cases 
with bright solitons: 

1. Hirota's equation (Hnls) @: ft : ft : ft = 1 : 6 : 0, 

qt + q X xx + G\q\ 2 q x = 0, (7) 

2. Sasa-Satsuma equation (SSnls) ft ' ft ■ ft = 1 : 6 : 3, 

<z* + fe^ + 6\q\ 2 q x + 3q\q 2 \x = 0. (8) 

Here the scaling convention mentioned above has been assumed and the eti terms eliminated. 

Some non-integrable special cases of Q have also been studied in the literature, including [jJE3 : ft : ft : ft = 1 : 
6:6, 

qt + q xxx + a{q\q 2 \)x — 0. (9) 



C. Multicomponent generalizations 

Both Hnls and SSnls allow various kinds of multicomponent generalizations, some of them integrable. The results 
of a Painleve analysis [Tlj can be summarized as follows: 

Casel: f qu + *~ + + = J (10) 

\ 92t + fexs + 3(|gi| 2 + l^nfec = 0, v 7 

which can be interpreted as a real 4-component modified Korteweg-de Vries (mKdV) equuation, reducing to Hnls for 
gi = q 2 etc., 

Case 2 Ea : ( + + ^j^^ + 1^!^^ + ? £?1 ^ 2 |2^ = 2' (11) 



and 



r« a p-?m. / 9it + 9i*M + 3(ki| 2 + ld 2 )fe + §<?i(|<?i| 2 + Id 2 )* = 0, . . 

Uase3[iJ. I ?2t + + 3( | ft)a + |, 2 | 2)52x + |, 2(kl |2 + {qA = 0j l^i 
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which reduce to SSnls under the above reduction, and a mixed case 

a 171172 . / ^i* + + f (|gi| 2 + |g 2 | 2 )gi^ + f giCQr^i^ + ^5^2^) = 0, , . 

Case4[MLU- { q2t + q2xxx +Z(\ qi \2 + \ q2 \*) q2x + Zq 2 ( q * qix + q * q2x ) = 0, {13 > 

which reduces to Hnls under q = qi = q 2 , a = 3 and to SSnls under q = q\ = q 2l a = 6. In each case we must of 
course adjoin the complex conjugated equations. Cases 1-3 are invariant under q 2 <-* q 2l while case 4 changes to the 
alternative form 

Case 4'- { Qlt + qixxx + % + \ q2 ^ qix + I + <&Q%x) = °> (u] 

\ Q2t + Qixxx + f (kil 2 + \q2?)q2x + ^(qiqL + q* 2 q2x) = 0. y > 

Under the reduction q 2 = cases 1,2,4 reduce to Hnls, and 3 to SSnls. 



D. The mKdV limit 



With complex and multicomponent equations it is important to make the following observation: we can always 
make the real, one-component reduction 

<ji(x,t) = Ciii(x,t),Vi, (15) 

where it is a real function and Cj are arbitrary complex constants. As a result of this all the equations mentioned 
before (and many others, including non-integrable ones) reduce to the real mKdV equation 

U t = u xxx + ku 2 u x . (16) 

(Note that for case 2 we need |ci| = \c 2 \.) This was observed already in 0, see eqs (21-25). A consequence of this 
rather simple observation is the following: 

All these complex and/or multicomponent systems always have multi-soliton solutions of the real mKdV type, with 
frozen complex and/ or multicomponent freedom. 

In the usual real one-component case the existence of multi-soliton solutions is a signal of integrability, but from the 
above we can see that this is not necessarily true in the complex/multicomponent case. In general it is essential that 
the individual solitons, from which the multi-soliton solution is built, are each allowed to have their own freedom 
of initial position and overall phase. That is, even if a one-soliton solution is of the type l|15fl . in the multi-soliton 
case each component soliton must be allowed to have its own parameters, including the complex coefficient(s) c;. 
Furthermore, under during scattering some of these parameters can change [Tlj . 

Thus in practice the reduction (|15fl trivializes the equation and the resulting solutions are hardly of interest. 
Nevertheless it seems that several recent studies have fallen into this trap and produced no solutions with genuine 
multicomponent or complex structure. This is quite evident from the proposed final results: for example the multi- 
component structure is trivialized into a constant factor in (see (3.15-16) or (3.20) or (3.25-26,32-33)), 01 ( see 
(17) or (24) or (27)) and ( see (10) or (13)), whereas the solutions are obviously real (after the gauge has been 
fixed) in Q (see Sec. Ill), ( see (2,3,15,16)) and (in Sec IV ki,r]i are real and H/G a constant). Below we will 
show that the reason for this often lies in the incorrect bilinearization that was used. 



E. Traveling-wave solutions 

Let us now return to the one-component equations (|7I8|I and consider their one-soliton solutions. For the purpose 
of orientation, let us first consider Hnls J7J). The usual traveling- wave ansatz 

q(x, t) = e Hax+U+u,)^ x +dt + S ), (17) 

where / is a real function (soliton envelope), leads to a pair of real equations, which are compatible, if 

b = a(3d-8a 2 ), (18) 

and in that case the solution can be parameterized as follows 

q (x, t) = eW^-Sc 2 )*^)^ c 

v ' cosh[c(a; + (3a 2 - c 2 )t + 6)] y ' 
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We observe that there are four free real parameters: a and c, which relate to the size and velocity of the soliton, and 
u),8 which give the constant complex phase and soliton position, respectively. 

If we use the same ansatz 1(17(1 in (JjJ we find that it works only under the additional condition a = 0, leading to 



q{x,t) 



\/2cosh[c(x - cH + 5)} ' 



(20) 



Since one parameter was lost the solution (|20|l is not general enough. Indeed, Sasa and Satsuma have derived a 
complex traveling- wave solution to , which does not fit to the usual ansatz 117fl , but has the form 



q(x,t) 



ia(x-\-(a —3c )t-\-oj) 



2eMe 2 '' + «Q 
e 4l > + 2e 2r > + Id 



a + ic 



, rj = c(x + (3a 2 — c 2 )t + 5). 



(21) 



We note that this has similar x, i-dependence as 1(19(1 but the functional form is different, also in the limit o->0, i.e., 
k — ► we obtain the real limit 1)20(1. 

It turns out that (|21|) is still not the most general one-soliton solution for this system, it is given by q = G/F, where 



G 
F 



rye' 1 



P e 

UI2 



\2p 2 + 2p* 2 



P_ e >7+V 



(P+P*) 2 



2p 2 



2p* 



4|p| 4 



2(77+r;*) 



= 1 



7e 1 



/re'* 



m = (h\ 2 P-\ P \Y) 



P 



Hh\ 2 + \p\ 2 ) 



(p- 




2 


me 2ri 


(p- 




+ 


2p 2 



px — p t + rj 



(0) 



(P + P*) 2 ' 
P, P, 7 



and r;^ ^ complex. 



(22) 
(23) 

(24) 
(25) 



Comparing with the original 1-soliton solution 1(21(1 . we have two extra parameters 7 and p. By 77 translation one 
finds that only p/7 matters, and if it vanishes we have the usual SS-solution, so this is a genuine new parameter. 
This parameter controls the oscillation, which appears not only in the carrier but also in the envelope (but in any 
case F > 1 so the solution is not singular). This solution was already obtained by Mihalache et. al. |2(j using inverse 
scattering transform, below we will derive it using the bilinear method. It is not easy to derive such a solution from 
a (complex) traveling wave ansatz, and Hirota's bilinear method is easier to use in this case. 



F. Outline 



In this paper we first present in Sec. [H] the bilinearizations that work and the one-soliton solution that is obtained 
by the expansion method. The detailed derivations and multi-solution solutions are made in Sec. IIIII 

It is well known that soliton equations can be organized into infinite hierarchies as described by the Sato theory pl| 
and that particular equations can be obtained from these hierarchies by various reductions. Indeed, one cannot have 
a full understanding of an integrable equation before its relation to integrable hierarchies is described. In Section ITTT1 
we will give the full picture by showing that the Sasa-Satsuma equation can be obtained from the general Sato theory 
as a reduction of the three-component KP hierarchy. The reduction can be made in two different ways producing 
two different bilinearizations. As intermediate steps of the reduction process we get either a (2 + l)-dimensional 
generalization or a complex 2-component generalization of the Sasa-Satsuma equation. 



II. DIRECT BILINEARIZATION AND ONE-SOLITON SOLUTIONS 



One can attempt to bilinearize the generic equation 

qt + q X xx + 6\q\ 2 q x + fiq{\q\ 2 )x = 0, 

with the standard substitution 



q 



G 
1' 



(26) 



(27) 
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where F is taken to be real and G complex. This leads to the equation 

F 2 [(D 3 X + D t )G-F] - f3GF(D x G-G*) - 3{D X G-F)[D 2 X F-F - |(/3 + 3)|G| 2 ] = 0, (28) 

which is quartic in F, G. We can see that if (3 = (which is the Hnls case) the equation splits naturally into two 
bilinear ones, (D x + D t )G-F = 0, and D 2 F-F — 2|G| 2 . In the general case (that includes the SSnls equation at j3 = 3) 
we could take 



DlF-F=Mp + 3)\Gf 



(29) 



as one of the equations, which leaves a trilinear equation 



F[{D Z X + D t )G-F] - (3G[D X G-G*\ = 0. 



(30) 



One might be tempted to require that in 130fl the terms in square brackets vanish separately, but this is not correct 
(as was already noted in £|), because it would result in more independent equations than there are unknowns and 
in effect force reduction to the real mKdV. [Clearly D X G-G* = d x (G/G*) — and therefore the phase of G has 
no x-dependence, and when G = R(x, t)e* 6 ^ is substituted into the remaining equation one finds that 0(t) must be 
constant, i.e., (|15() .] As a matter of fact, this sort of brute force bilinearization turns out to be precisely the reason for 
the trivialization of the complex/multicomponcnt freedom mentioned before. Unfortunately this incorrect approach 
has been used quite frequently, see, e^., [3 (3.6), P3 (9,29,44,47), (12), (26), Q (21-24), [l| (7c), 
(20), H| (43), (18), [13 (9,16), H3 (38). 

The trilinear equation 1|30|) can only be split into two bilinear ones by introducing a new dependent variable. There 
are two acceptable ways to do this, resulting in 



((6 
((6- 



1-209 
209 H 



D 2 X F-F 
f 3)D t )G-F 
3)D t )G*-F 
D 2 X G-F 
D 2 G*F- 



§G0 + 3)|G| J 

3{3D X H-F, 

3f3D x H*F, 

-HF, 

—H*F. 



(31) 



or 



D 2 X F-F 
(D x + D t )G-F 
(D x + D t )G*-F 
D T G-G* 



f(/? + 3)|G| J 
0SG, 
-(3SG*, 
SF, 



(32) 



where the new dependent variable has been called H and S, respectively. Note that S is pure imaginary, H complex, 
and that HG* — H*G = D X F-S. These splittings are acceptable, because they introduce equal number of new 
functions and new equations, and furthermore for integrable equations and soliton solutions the new functions turn 
out to be expressible in terms of polynomials of exponentials. Thus for any (3 we can give for l|26|) a bilinear form in 
terms of three bilinear equations for three functions, but it should be emphasized that the fact that an equation can 
be written in a bilinear form does not by itself imply that the equation is integrable, or that the new functions S, H 
are r- functions, although it is the case when (3 = 3. 

The one-soliton solution can be obtained as usual by substituting the expansions 



1 + e 2 F 2 + e 4 F 4 



G = eGi + e 3 G 3 



(33) 



accompanied with suitable ansatze H or S, into (|31() or l|32() . and truncating at some power of the formal expansion 
parameter e. For Hnls {(3 — 0) the expansion can be truncated by keeping terms up to e 2 , but for SSnls (/3 = 3) we 
must go up to e 4 obtaining F,G as given in I|22I25|I . with the auxiliary functions 



S 



(p-p*)(l7l 2 -MV +rr , 



H = —"ip e ' — p p e ' 



m 



*2 * 2 

_fP_ e 2 n +rf + P_P_ e r ) +2r ) * 
p2 p*2 



(34) 
(35) 



which also are polynomials of exponentials. It is not known whether the expansion can be truncated for any other 
value of f3. 
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III. THE SASA-SATSUMA EQUATION AS A REDUCTION OF THE 3-COMPONENT KP 

HIERARCHY 

We will next explain how the Sasa-Satsuma equation and its multi-soliton solutions can be obtained from the 3- 
component KP hierarchy by suitable reductions. It turns out that there are two different reduction routes leading to 
the Sasa-Satsuma equation; both are two-step reductions but the intermediate equations are different. We will first 
describe the starting point (3-component KP hierarchy) and then the two kinds of reductions. 



A. 3-component KP hierarchy and its r-functions 

In general, the 3-component KP hierarchy has r-functions depending on three infinite sets of variables x = 
x, X2, X3, . . . , y = y, 1)2, 2/3, ■ ■ ■ and z = z, Z2, Z3, . . . , and is defined in terms of vector "eigenfunctions" </>(x), i/>(y) 
and x(z) and "adjoint eigenfunctions" </>(x), ^(y) and x(z). We should emphasize that, at this point, these six 
eigenfunctions are independent of one another. In general, they are only assumed to satisfy the linear equations (for 
n = 2,3,...) 



&„0 = ^0, -d Xn 4> = {-d x ) n fr 

dyji = d^, -dyj> = (-By)"*, 



(36) 
(37) 
(38) 



Here we consider the special case in which only dependence on x, X2, X3, y and z is active and so the vectors <f>(x, X2, £3) 
and 4>{x, X2, X3) satisfy 



d X2 <t> = dl<j>, d X3 (j) = dl<j), d X2 <p = -d%.<t>, d X3 <j> = dl<j>, 

and ip(y), x( z )i ^(y) an d x( z ) are arbitrary vector functions of a single variable. 
A potential matrix m is defined by 

d x m = (j><t>, d y m = ipxp* , d z m = xx\ 

which can be integrated to 

m = c + J 4><fi t dx + J tp'i(} t dy + J X^dz, 

where c is a constant matrix. As a consequence of l|36|) . we also have 

d X2 m = 4> x 4> - <f>4>l, d X3 m = § X J> - §J> X + <f>4> xx , 



Now define the r-functions 



/ = M 



(39) 

(40) 
(41) 

(42) 
(43) 





m 







m 





9 = 







. 9 = 








and 



h = 



m <fi 

-x l 



h = 



m x 
-4> 



(44) 



(45) 
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{D 2 x ~D X2 )g 


f 


= o, 


{D 2 x +D X2 )g-f = 0, 










{D 2 x -D X2 )h 


f 


= 0, 


(D 2 x +D X2 )h-f = 0, 








(Dl- 


YiD x D X2 -AD x ,)g 


f 


= o, 


{Dl ~ 3D X D X2 - AD, 




/ 


= 0, 




-3D x D X2 -4D x Jh 


f 


= 0, 


{Dl - 3D X D X2 - AD, 




/ 


= o, 




DyD x f ■ f = 






D z D x f ■ f = -2hh, 










D y (D 2 x ~D X2 )g 


f 


= 0, 


D y {D 2 x + D X2 )g ■ f = 


o, 








D z {Dl-D X2 )h 


f 


-o, 


D z (D 2 x + D X2 )h-f = 


0. 







By considering Jacobi determinantal identities involving /, g, g, h and h and their derivatives with respect to x, x%, 
x$, y and z, one may compile a complete list of bilinear equations that are satisfied by these functions. The bilinear 
equations given below are the only ones from this list that will actually be used in the rest of this paper: 

(46) 
(47) 
(48) 
(49) 
(50) 
(51) 
(52) 

As is typical for the multicomponent KP hierarchy, this set of equations appears to be over-determined as it stands, 
having many more equations than dependent variables. But we already know that it has a rather general set of 
solutions given above (even containing several arbitrary functions of one variable). It turns out that there exist 
exactly the right number of differential relations amongst these equations to guarantee their compatibility. There 
is some freedom in choosing the primary or independent equations, one choice is ()46l47l48fa , l50|) [7 equations for 5 
functions and two dummy independent variables] . The remaining equations are consequences from these or possibly 
just restrict some constants of integration. As an example consider (|48b ). From (|46b l we can determine g X2 , from 
BHp) .9x3 an d from their cross derivatives we get an equation for /. But now doing the same computation for g from 
lj%ob ) and ljl£h) we get the same equation for / and thus J3Bb) does not add essential information. 
In order to write (|46I52I) in nonlinear form let us first introduce the dependent variables 

?=~, q = j, r = j, f=~ and $ = | (log J% , * = i(log/) X2 . (53) 

Converting the bilinear equations into nonlinear form and eliminating dependence on the auxiliary variable X2 one 
obtains 

Qxxx + 6q x $x + 3q{®xx + - Qx 3 = 0, 

(jxxx + 6q~x$x + 3q{®xx - - q X3 = 0, ^ 
r X xx + 6r x $ x + 3r{<5> xx + * x ) - r X3 = 0, 1 ' 



together with 



r xxx + 6r x $ x + 3r($ xx - * x ) - r X3 = 0, 



® v = -\qq, $ z = -hrr, 



*y = ~h(QxQ- qqx), * z = -\{r x r - rr x ). 



Although this looks superficially like a (3+l)-dimensional system, it in fact describes a family of (2+l)-dimensional 
systems. The y- and z-dependence arises in such a way that it could be replaced by single variable corresponding to 
any linear combination of y and z. 

In the next sections we will describe a two stage reduction of this system in which a calculation similar to that used 
above will give the Sasa-Satsuma equation. 

B. First reduction, step 1 

The previous set of equations contains two dummy variables, x-i and one of y, z. In this reduction we will eliminate 
the dummy variables by keeping just the leading terms in x-i and y — z. 

We start by considering eigenfunctions and adjoint eigenfunctions possessing the symmetry 

4>{x,-x 2 ,x 3 ) = 4>{x,x 2 ,x 3 ), (56) 

and the other eigenfunctions having pairwise identical forms: 

tp{a)=x{a), x(a)=i>(a)- (57) 
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This reduction may be shown to be a natural generalization of the 3-component version of the C-reduction described 
in H2. 

Now we explore the consequences of this symmetry on the r-functions. Using the independent variables y = £ + rj 
and z = £ — 77, the symmetry (|57H gives 

0(y) = xi$ + v) =x(z)\ n ^_ n , x(z) = ">!>(£ -v) = ^(y)\ n ^- n - 

For the potential to it is then easy to see that 

m(x, -x 2 ,x 3 , -77) = m t (x, x 2 ,x 3 ,£„ <rj), 
as long as the constant matrix c in (|41fl is taken to be symmetric. Hence 

f(x, -X2,X 3 ,£, -7]) = f(x, X 2 ,X 3 ,£, 7?). 

In a similar way, 

g(x, -x 2 ,x 3 , £, -77) = h(x, x 2 ,x 3 ,£, 77), h(x, -x 2 ,x 3 ,£, -77) = g(x, x 2 ,x 3 ,£, 77). 
Next we consider the Taylor expansions of the eigenfunctions with respect to x 2 and 77 and obtain 

<f>(x, x 2 ,x 3 ) = (f>(x, 0, x 3 ) + x 2 cf) xx (x, 0, x 3 ) + 0(x\), 

while the symmetry (|56|l gives 

4>(x, x 2 ,x 3 ) = 0(x, 0, x 3 ) ~ x 2 cj) xx (x, 0, x 3 ) + 0{x\). 

By a similar argument 

^(y) = rP(0 + o(v), x(z) = x(0 + O(v), 0(y) = x(0 + 0(77), = 0(0 + 0(77). 

For the potential to the expansion is 

m(x,x 2 ,x 3 ,£,ii) = m(x,Q,x 3 ,£,0) + x 2 (0J.0* - 00^) (a;, 0,2:3) + 0{x\,r]). 
For the r-functions we then get 



/ = /o(»,»3,0 + 0(xl,f]x 2 ,T] 2 ), 



and 



/o = |mo|, 

where too = m(x, 0, x 3 , £, 0) satisfies 

mo.x = (t> 4>l, rn , X3 = <t>o,xx<t>o ~ 0o,s0o,; C + 0000,* > ™o,£ = 0X* + X0*, 
and O — 0( ;e j 0; ^3) satisfies the single linear pde 

9*3 00 = 9 *0o- 

Also, 



where 



5o 



0(xl,i]) 


ft = 


7 + x 2 h 2 


to 0o 




to o 


-x* 


, fto = 


-0* 



and 



.92 



m 00,** 

-x* 



m 00 00.* 

-0o 
-X 4 



too Oj2; 
-0* 



TOo 0o 0o. ; 

-0o 
-0* 



(58) 
(59) 
(60) 
(61) 
(62) 
(63) 
(64) 
(65) 
(66) 
(67) 
(68) 
(69) 

(70) 
(71) 
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Finally, because of 160|) 

g = h - x 2 h 2 + 0{xl,rj), h = g Q - x 2 g 2 + 0{x\,rf). (72) 

The above discussion shows that, up to leading orders in x 2 and 77, the original five r-functions /, g, h, g, h can be 
written in terms of the five r-functions fo, 30, 32, ho, h 2 depending only on x, X3 and £. 

The final part of calculation is to identify an appropriate set of five bilinear equations involving these five r-functions. 
Applying the reduction to l|46[) -H50 [) gives 

D l9o ■ fo = 32/0, D 2 x ho ■ fo = h 2 fo, (73) 
{Dl - 4D X3 )g ■ fo = -3D x g 2 ■ f , (D x - 4D X3 )h ■ f = -3D x h 2 ■ /„, (74) 
D c D x f ■ fo = -4g h . (75) 

Notice that in this reduction, (|4T1 and (|4TI) both give J73J), |@HJ and |@j5J| both give l(74*jl and the sum of the equations 
in gives l[7H)l. 

This set of bilinear equations (|75)) - lf75'|) is the Hirota form of a (2 + l)-dimensional Sasa-Satsuma equation. If we 
define 



9° ^° 52 h 2 1 . 

?=T> r = T' Q2 = -r, r 2 = — and $ = 5 (log Jo)* (76) 
Jo Jo Jo Jo 



then lUJJ give 
JH} give 
and |j75)| gives 



<72 = + 4q$ x , r 2 = r xx + 4r$ x , (77) 
+ 12g x $ x - 4^3 = -3g 2 ,x, r l;cl + 12r x $ x - Ar X3 = -3r 2 , x - (78) 
$ c = -gr. (79) 



After eliminating q 2 and r 2 one obtains 



Gqx^x + iq^xx = q X3 , 

6r x $ x + 3r$ xx = r X3 , (80) 
= -gr. 



If we now set r = — q*, X3 = —t and use U = Q x we get a (2 + l)-dimensional Sasa-Satsuma equation 

qt + q X xx + 6q x U + 3qU x = 0, 

U( = (\q\ 2 )x- 



(81) 



C. First reduction, step 2 

In order to make a dimensional reduction from (2 + 1)- to (1 + l)-dimensional, we make a further rotation of 
coordinates 

x = i(X + S), £ = i(X-H), (82) 

and then choose eigenfunctions so that the r-functions will be independent of S. Then both X and H derivatives in 
(|SU|) become x derivatives and we obtain the Sasa-Satsuma equation with two complex fields 

qxxx - &q x qr - 3q(qr) x - q X3 = 0, r XX x - 6r x qr - 3r(qr) x - r X3 = 0. (83) 

In order to keep the solution structure in this dimensional reduction, it is necessary to choose eigenfunctions </>, ■»/> 
and x so that they are separable with a common dependence on S. A natural way to achieve this is to take 

(0„)< = \ie PiX+p * x * -> Keh-eh^+v^ , (84) 
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^= foe-Pit ^ Ai e2^ s e -2^, 



where Xi, fii, i>i and are constants. As a result of this choice of eigenfunctions we have 

1 , _ x„ , ,„3 , _3x_ i, 



(85) 
(86) 



(m )y = Cij + e2 Pi 



Pt + Pj 



i „ 



iW a , hvpi+ Pj ^0, (87) 



where cy are constants of integration [if pi +pj =0 we must choose the coefficients properly so that from (|41|) we get 
a constant, which can then be absorbed into the c-matrix]. As a consequence of the C-reduction, the matrix cy has 
to be symmetric. In order that /o = |toq| be independent of S we must have 



e2 Pi " = constant 



(88) 



and, for each i,j G {1, . . . , L}, 



dje 2 



(pi+p j)E 



constant. 



(89) 



These are satisfied if and only if Yli=iPi = an d f° r each i,j g {1,...,L} either cy =0 or ^ + Pj = 0. 
Consequently, we take L = 2N and then 



PN+i = —Pi, Vi = 1, . . .,N, dj = Si+Nj Ci - Sij+NCj, Vi, j £ {1, . . . ,2N}. 



(90) 



Finally, we show how to obtain solutions of the usual Sasa-Satsuma equation JSJ, in which r = — q* , where * stands 
for complex conjugation. In order that this comes about, we must have f Q real and = —go, h* 2 = —g2, and so we 
must impose the relations 



<i>* = p<t> , r = Px, x* = p*j>, 



(91) 



where P is a permutation matrix. The simplest realization of these conditions is to take N = 2M, choose the 
permutation to be 





E 








I 




















E 


Vo 





I 


o 



with M x M blocks, and 



</>o = 



\ M e2 PMX+p ™ X3 

\* e 2Pi X+ Pi 3x3 



J- * * 3 

\* p2Pm x +Pm x 3 



( o \ 








Mie2 



PiX 



i , v 

* 77 p, x 



\v* M e?^ x ) 



I o \ 






and x 



(; 



PiX 



VMS' 1 

tie? 



PmX 

pIX 



(92) 
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where we have changed the notation for coefficients in order to conform with (|91|l . Since the e 2 Pk "-factors in (|84I87|) 
will eventually cancel out with the above choices we do not include them in these formulae, but in order to compensate 
this omission we must replace <9™</>o by (2dx) n 4>o, e.g., in J7TJ. 

Taking all constants of integration c$ = 1 in (|90|) gives the M-soliton solution. In particular the one-soliton solution 
shown in I|22I25[I is obtained for M = 1, Ai = 1,/ii = — p, v\ = —7. 

If we set x 3 = —T and following (|%2"|l replace D x and with Dx in the bilinear equations (|73I75|) they become 



D 2 x 9o ■ fa = 52/0, 
(D% + AD T )g ■ f = -3D x g 2 ■ fo, 



D 2 xh • fo 



(93) 
(94) 
(95) 



This is the same as Ij.HII for /3 = 3, if we identify /o = F, go — ~G and g 2 = H. The multisoliton solutions are obtained 
from (|87I66I70I71() with 192(1 . [But please remember that due to the simplified expressions we have d x <fio — (2dx) n 4>o-] 

D. Second reduction, step 1 

To obtain the alternative bilinear form of SSnls we carry out the reduction process in a different manner. This 
process will take us via a "coupled Sasa-Satsuma equation" as opposed to the (2 + l)-dimensional Sasa-Satsuma 
equation. 

Firstly we need to introduce two new r-functions, 





m 





<f>x 




m 




X 


s = 










, s = 


-~i 




















~<t>x 









(96) 



In addition to I|4t)l52|l we now have some further bilinear equations satisfied by these r-functions together with the 
original five r-functions (|43I45() : 



D Z (D 2 X ~ D X2 )g ■ f = Ask 
D y (D 2 x -D X2 )h-f = -Asg, 
D x h ■ g = sf, 



D z (D 2 x + D X2 )g- f = 4sh, 
D y {D 2 x + D X2 )h-f = -Asg, 
D x h- g — sf. 



(97) 
(98) 
(99) 



Again these equations are not all independent. Altogether there are 7 dependent variables f,g,g,h,h,s,s and two 
dummy independent variables and therefore we need 9 independent equations. We can take, e.g., 1)461 1471 l4*8b . 1501 199JI . 
and then the other equations are consequences of these. [In practice it is best to keep the full set at one's disposal.] 
If we change variables to £ and r\ using 



y = £ + ri, 



(100) 



then taking sums and differences of some of the equations, for instance 1151b ,) and (|97b ) we get some equations 
containing only ^-derivatives and others containg ^-derivatives. In the following we will only use the ones containing 
^-derivatives, they are 



Dz (Pi - D X2 )g ■ f = Ash, D^Dl+D x ,)g-f = Ash, 



D S (D 2 X -D X2 )h-f = -Asg, D^(D 2 X + D X2 )h ■ f 
D £ D x f-f = -2(gg + hh). 



This leaves us with equations H46l49ll9^1H01ll03|) . 

At this stage we have not yet carried out a reduction. If we now do a second change of variables 

x=i(X + E), £=i(X-S), 



(101) 
(102) 
(103) 



(104) 



we can achieve a dimensional reduction in a manner similar to the dimensional reduction in Section IlII CI i.e., by 
expanding in 77, H and keeping only the leading terms. After also eliminating the x 2 dependence, and denoting 
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X3 = —T, we finally obtain the following set of equations 

(D% + D T )g ■ f = 3sh, (D x + D T )g ■ f = 3sh, (105) 

(D x +D T )h-f = -3sg, (D x +D T )h-f = -3sg, (106) 

D x h-g = sf, D x h-g = sf, (107) 

D 2 x f-f = -2(gg + hh). (108) 

This is a coupled Sasa-Satsuma equation with complex fields. The nonlinear form obtained with the substitutions 

9-9 h _ h 

" 1 " 1 ' 1 ' J 

is 

qr + qxxx - Gqxqq - 3r(qr) x = , 

qr + qxxx - Qqxqq - 3r(qf) x = , , log > 

tt + rxxx — Qrxrf — 3q(rq)x = , 
?t + rxxx - 6f x rr - 3q{fq) x = , 

which was proposed already in 01 ■ H we take q — —q*, f — —r* we obtain (|13|) . 



E. Second reduction, step 2 



The final reduction on this system is a reduction of C type, this is obtained as in the other bilinearization by 
identifying 

g = h, h = g, s= -s. (110) 

This gives us the alternative bilinear form of the Sasa-Satsuma equations: 

(D x +D T )g-f = 3sg, (D x + D T )g ■ f = -3sg, (111) 
Dxg ■ 9 - sf, (112) 

D 2 x f ■ f = -4gg. (113) 

With f — F, g — — G, g — G* , s — S these equations yield (j32(l for [3 = 3. The solutions for this alternative form will 
be the same as in the first case and the nonlinear form of these equations is (j83|l . with r replaced by q. This bilinear 
form of the system requires a single pure imaginary auxiliary variable s, whilst the other bilinearization involves a 
complex auxiliary field h, and consequently here we need four bilinear equations rather than five. 



IV. CONCLUSIONS 



In this paper we have shown how the Sasa-Satsuma equation fits into the general theory as a reduction of the 
three-component KP hierarchy. As a result we have obtained two possible bilinearization for SSnls, 1)31(1 and l|32Jl . 
and formulae for constructing multisoliton solutions, (|87l 1921 Ififil 1701 fTTI I96|l . In the reduction process we have also 
obtained two intermediate equations, JHTJ, and H109|) . of which the (2 + l)-dimensional (|81|l seems to be new. 
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